Kinks in the electronic dispersion of the Hubbard model away from half filling 
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We study kinks in the electronic dispersion of a generic strongly correlated system by dynamic 
mean-field theory (DMFT). The focus is on doped systems away from particle- hole symmetry where 
valence fluctuations matter potentially. Three different algorithms are compared to asses their 
strengths and weaknesses as well as to clearly distinguish physical features from algorithmic artefacts. 
Our findings extend a view previously established for half-filled systems where kinks refiect the 
coupling of the fermionic quasiparticles to emergent collective modes, which are identified here as 
spin fluctuations. Kinks are observed when strong spin fluctuations are present and, additionally, a 
separation of energy scales for spin and charge excitations exists. Both criteria are met by strongly 
correlated systems close to a Mott-insulator transition. The energies of the kinks and their doping 
dependence fit well to the kinks in the cuprates, which is surprising in view of the spatial correlations 
neglected by DMFT. 

PACS numbers: 71. 27.-|-a,71.30.+h,74.25.Jb, 71.28. +d 



I. INTRODUCTION 

Collective bosonic modes can modify the low-energy 
electronic properties on the characteristic energy scale of 
these bosonic excitations. Well understood are the signif- 
icant mass enhancemenlPJ and the kinks in the electronic 
dispersiorP in materials with strong electron-phonon cou- 
pling. These kinks are abrupt changes of the slope of the 
electronic dispersion which occur at energies of the or- 
der of the Debye energy which is much smaller than the 
Fermi energy in such systems. 

In this paper we study kinks in the electronic dis- 
persion which are related to emergent collective bosonic 
modes of the purely electronic system rather than by 
a coupling to external bosons such as phonons. It 
was shown by Byczuk et al. in the framework of dy- 
namic mean-field theory^ISl (DMFT) that such kinks are 
a generic feature of strongly interacting systems and re- 
quire no externally coupled bosons.^' Their argument re- 
sides on the structure of the DMFT equations and the 
spectral density of the local propagator. 

Our work supplements the mathematical argument in 
Ref. [5] by the physical picture that the internally gener- 
ated emergent collective modes provide the bosons which 
generate the kinks in the dispersion. This point of view 
has been already put forward for the half-filled case where 
the physical situation is particularly clear.El We will show 
that in the model studied below these collective modes 
are spin fluctuations of local moments formed due to the 
strong local Coulomb interaction. In the Mott insulating 
phase, stable local moments are generated which tend to 
form magnetically ordered phases due to residual spin- 
spin interactions. In the metallic phase, however, these 
local moments are completely screened at zero energy and 
the spin fluctuations at intermediate energies live only for 
a finite time because they decay into particle-hole pairs 
(Landau damping) . We call the spin fluctuations appear- 
ing as resonances emergent modes because they become 
long-lived for even larger interaction in the insulating 



phase. Photoemission experiments support our view in 
antiferromagneticP'21 and in ferromagnetic systems.^^ 

As will be shown, there exist two conditions for the 
kinks to appear in the dispersion relation. For one, strong 
emergent collective modes (spin fluctuations) must be 
present in the system. This is signaled by a pronounced 
maximum in the imaginary part of the spin suscepti- 
bility at a finite energy Wmax- For another, the en- 
ergy scale for charge excitations Wchargo must be much 
larger than the characteristic spin-excitation energy, i.e. 
|i^chaigc| > l^inaxl- Here, the energy scale for charge 
excitations is set either by a single-particle energy or 
the characteristic scale for collective charge excitations, 
depending on which is smaller. Both criteria are met 
in strongly correlated metallic systems close to a Mott- 
insulator transition. 

The term "collective bosonic modes" in solid-state sys- 
tems usually refers to dispersive, i.e., wave- vector de- 
pendent, excitations such as (para-)magnons or charge 
density waves. Such modes are characterized by the 
momentum dependent spin- and charge-susceptibilities 
Xmag(g,w) and Xchargclg, w), respectively. It may be 
put forward that in DMFT no true collective modes 
exist because they do not propagate properly. In- 
deed, in the strict limit of infinite dimensions collec- 
tive two-particle properties are local because they are 
momentum-independent except for some momenta of 
measure zero.l^i^illYet we do not adopt this puristic view 
and stick to the wording of collective bosonic modes for 
three reasons. First, spin- and charge fluctuations are 
two-particle properties and as such comprise at least two 
elementary fermionic excitations so that they are collec- 
tive and bosonic in this sense. Second, even in the strict 
DMFT, the spin- and charge-susceptibilities are true lat- 
tice quantities involving all lattice sites. Otherwise, their 
imaginary parts would show only one or two S peaks in- 
stead of broad continua, see for instance Ref. [151 Third, 
generically the DMFT is taken as an approximation to 
finite dimensional systems. Then, the susceptibilities ac- 



quire a non-trivial g-dependence through the inner mo- 
mentum and frequency sums in the Bethe-Salpeter equa- 
tions, even though t he irr educible two-particle vertex is 
assumed to be localP^^^^. In the present work, we ex- 
tract the characteristic energy scales relevant for the dis- 
persive collective modes from the local susceptibilities 
x(a;) = jfJ2qX{Q':^) because they readily reflect these 
energy scales and can be obtained in DMFT from the 
effective impurity model. 



We investigate the low-temperature phase of the Hub- 
bard model away from half filling using the DMFT with 
three different impurity solvers in order to be able to 
clearly separate physical features from possible algorith- 
mic artifacts. We employ two numerical renormaliza- 
tion schemes, namely the density-matrix rcnormalization 
group-^' (DMRG) and numerical rcnormalization groupl- 
(NRG) as well as a diagrammatic approach, th e so-ca lled 
enhanced non-crossing approximation (ENCA) .EflEll xhe 
results of all three algorithms agree very well. The re- 
maining differences can be understood on the basis of 
the known strengths and weaknesses of the approaches. 



The prevailing effect of the interaction is a signifi- 
cant renormalization of the bare electron mass due to 
correlations. The Fermi liquid theory^H connects this 
mass renormalization to a fundamental energy scale T* 
below which the Fermi liquid picture of renormalized 
weakly interacting, long-lived quasiparticles applies. In 
strongly correlated electron systems the properties of 
these renormalized quasiparticles may differ dramatically 
from those of non-inte ractin g electrons as exemplified in 
heavy fcrmion systems.^ ^^ * ^^ * 



The article is laid out as follows. The model and the 
relevant theoretical background are described Sec. |lT] In 
particular, the basic ideas of all three impurity solvers 
are discussed. In order to gauge their differences a com- 
parison of the single-particle dynamics is shown in Sec. 
|III| We begin with the results for the self-energies. Then 
we turn to the extraction of the Fermi liquid low energy 
scale T* and relate it to the maxima in the local charge 
and spin susceptibilities, in Sec. |IIIB[ We will explicitly 
show that T* also occurs as energy scale in the dynamic 
spin susceptibility but not in the local charge suscepti- 
bility. Therefore, the connection made between emer- 
gent spin-fluctuations and the the kinks in the electronic 
dispersiorP extends to the metallic regime away from half 
filling. This main results of our work is contained in Sec. 
|IV| where also the doping dependence of the kink energies 
is compared to experimental results. Finally, our findings 
are summarized in Sect. fVl 



II. MODEL AND METHODS 

A. Hubbard Model and Dynamic Mean-Field 
Theory 

We consider the single-band Hubbard model which is 
the simplest model for correlated electrons on a lattice 






^E 



ni-friii 



U 



E^ 



(1) 



M > n. 



The operators c]^ (q^) create (annihilate) electrons with 

spin a = {t, i} at lattice sites i and riicr = cl^c^^ is 
the occupation number operator. The first term in the 
Hamiltonian M describes electronic hopping with ampli- 
tude t, where {i,j) indicates nearest neighbor sites. The 
local Coulomb repulsion is incorporated in the second 
term with matrix element U. The last term sets the lo- 
cal single-particle levels and includes the chemical poten- 
tial in a way, that fi = corresponds to an electron-hole 
symmetric situation, i.e., half-filling. 

Despite its simplicity, the exact solution of the Hub- 
bard Hamiltonian (fll) has only been possible in one 
spatial dimension, ^^'^ for a recent book on this topic 
see Ref. 1271 In order to obtain an approximate solu- 
tion for higher dimensional syste ms w e employ the dy- 
namic mean-field theory (DMFT),EB] for recent reviews 
see Refs. [^50gi 

A non-trivial, but considerable, simplification is ob- 
tained in the limit of infinite coordination number (in- 
finite spatial dimension) if the hopping matrix elements 
are rescaled appropriatelyl^^HSIl jjj ^j^jg limit DMFT rep- 
resents the exact solution.lSBI When applied in finite di- 
mensions, the major approximation of DMFT consists in 
treating all non-local correlations in a mean-field man- 
ner while the correlated local dynamics is faithfully re- 
tained. This translates to the self-energy being momen- 
tum independent,'^^ S(fc, z) — >■ S(z) where we use z as 
a general complex energy argument with finite imaginary 
part. 

Then, the local Green function of the lattice problem 
reads 
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where N is the number of lattice sites and e.- the bare 

k 

dispersion. This local lattice propagator equals the local 
Green function of an effective single-impurity Anderson 
model (SIAM) 
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embedded in a medium characterized by the hybridiza- 
tion function V{z). Thus, the embedding medium is a 



dynamic medium and it is not independent from the so- 
lution of the SIAM, but it has to be determined self- 
consistently as in any generic mean-field approach. 

For a given guess for the hybridization function the lo- 
cal Green function is determined by a suitable numerical 
algorithm which is commonly referred to as the employed 
'impurity solver'. This yields the self-energy 
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(4) 



which in turn is used to obtain the local lattice propa- 
gator G{z) via Eq. ([2|. It is in this step that the lattice 
structure enters. The self-consistency cycle is closed by 
reorganizing Eq. m and to extract a new guess for T{z). 
This cycle is iterated until convergence is reached in the 
pragmatic sense of a tolerable deviation of two succes- 
sive results for the local Green function, the self-energy, 
or the hybridization function. 

The nontrivial step in this cycle is the solution of the 
effective SIAM. Due to the long history of DMFT, a 
multitude of different impurity solvers for treatin g th is 
model exists such as iterative perturbation theory,E2lEll 
exact diagonalizatiorPSESl or sev eral variations of quan- 
tum Monte-Carlo schemes !^^ l ^^" ^ In this work, we em- 
ploy the enhanced non-crossing approximation, the nu- 
merical renormalization group, and the dynamic density- 
matrix renormalization group as impurity solvers and 
compare their results. 



B. Enhanced non-crossing approximation 

The enha nced non-crossing approximation 
(ENCA),'22EiIli 42 sometimes also called one-cro ssing 
approximation, ^3 is a thermodynamically conserving^SSl 
approximation for the SIAM which utilizes the expansion 
with respect to the hybridization between the impurity 
electrons and the conduction band.l^^Ml It ext ends the 
usual non-crossing approximation (NCApfilSil to finite 
values of the Coulomb repulsion U via the incorporation 
of the lowest order vertex corrections, which are nec- 
essary to produce the correct Schrieffer- Wolff exchange 
coupling and the order of magnitude of the low energy 
Kondo scale Tk in the problem. The impurity spectral 
functiorPSl and dynamic susceptibilities^ are extracted 
directly for real frequencies without any adjustable 
parameters. 

The NCA is known to violate Fermi liquid proper- 
ties at very low temperatures and some pathological 
struc ture appears at the Fermi level below a pathology 
scale .E2E3] xhe ENCA removes the cusps associated with 
this pathology^*^ and significantly improves the Fermi 
liquid properti es of the spectral functions and dynamic 
susceptibilities .2322 However, the skeleton diagrams se- 
lected within the ENCA still suffer from an imbalance 
between charge and spin excitations. While magnetic 
properties, i.e., the magnetic susceptibility, are described 
excellently, charge fiuctuations are not as well accounted 



for.^^ In the spectral functions, where spin and charge 
fiuctuations contribute equally, this leads to an overes- 
timation of the height of the many-body resonance at 
the Fermi level at too low temperatures. The dynamic 
charge susceptibility is overestimated for very low fre- 
quencies. These shortcomings are related to the thresh- 
old exponents of the auxiliary ionic propagators and tend 
to be worse in parameter regimes with substantial valence 
fiuctuations, i.e., in the mixed valence regime or at small 
Coulomb repulsions U. 

Within the DMFT the overestimation of the many- 
body resonance of the impurity solution might lead to a 
violation of causality in the self-energy due to the sub- 
traction occurring in Eq. (HI) . In parameter regimes with 
considerable valence fluctuations, i.e., at small U or at 
large doping, this fact limits DMFT calculations to tem- 
peratures above the characteristic low temperature scale 
T* of the lattice. Because T* — )■ on approaching the 
Mott-Hubbard metal-to-insulator transition (MIT) the 
ENCA can be used as impurity solver down to very low 
temperatures in the vicinity of the MIT. 

Detailed comparisons of various approximations based 
on the hybridization expansion can be found in the liter- 
ature, e.g., in Refs. [51 HTl F551 and 1351 



C. Numerical renormalization group 

The NRG is a very powerful tool for accurately cal- 
culating equilibrium properties of complex quantum im- 
purities. Originally developed for treating the single- 
channel Kondo Hamiltonian,'^ this non-perturbative ap- 
proach was s ucces sfully extended to the Anderson im- 
purity modeP^ES and other more complex quantum 
impurities.'^ At the heart of this approach is a loga- 
rithmic discretization of the continuous conduction band, 
controlled by the discretization parameter^ A > 1. 
Using an appropriate unitary transformation, the dis- 
cretized Hamiltonian is mapped onto a semi-infinite 
chain, with the impurity coupled to the open end. The 
iVth link along the chain represents an exponentially de- 
creasing energy scale, ujn oc A~^'^. Using this hierarchy 
of scales, the sequence of finite-size Hamiltonians "Hat 
for the iV-site chain is solved iteratively, discarding the 
high-energy states at the end of each step to maintain 
a manageable number of states. The reduced basis set 
of Hn thus obtained is expected to faithfully describe 
the spectrum of the full Hamiltonian on a scale of wjv, 
corresponding to the temperature T/v ^ wat.!^ Because 
the thermal Boltzmann factors suppress the contribu- 
tions of higher lying energy states exponentially, the re- 
duced NRG basis set of Hn is sufficient for an accurate 
calculation of thermodynamic quantities at temperature 
Tn- 

Dynamical quantities, however, such as impurity 
Green functions and susceptibilities, require the infor- 
mation on all energy scales. In a recent exte nsion of the 
NRG to real-time dynamics out of equilibriunPSlSl] ^ com- 



plete basis set for a Wilson chain of length N has been 
identified which is a lso us ed for the accurate calculation 
of spectral functionsl^^MI Additionally, the discretization 
error is reduced by averaging over several different dis- 
cretizations of the conduction band.^** In order to re- 
duce the arbitrariness in the spectral broadeningSStSIlthe 
single-particle self-energy entering the DMFT equations 
is obtained from an exact expression of a ratio of two 
correlations functions!^ Since the local dynamic bosonic 
spin and charge susceptibilities are calculated directly 
from the NRG raw data^^ more pronounced broadening 
artifacts occur. 

In this work we use a discretization parameter A = 2 
and keep approximately 1700 states in each NRG itera- 
tion step. Eight different band discretizationgS^l are av- 

"651 



eraged and the artificial logarithmic broadening^^l varies 
between b = 0.08 and b = 0.2. 



for discrete complex frequencies ujj — ir/j, r]j > 0. In 
order to obtain G{uj) at the real axis with continuous 
spectral density p{uj) := ImG'(a;)/7r, we use the least-bias 
deconvolution algorithm. 

The key advantage of the correction vector DMRG is 
a good energy resolution for low and high frequencies 
uj. With correction vector DMRG combined with least- 
bias deconvolution, local Green functions^" ^^ and local 
susceptibilities^^ have been computed successfully. 

In this work, we use a fixed distance Aw = O.ID be- 
tween two successive frequencies Wj and w^+i, the arti- 
ficial broadening is rjj = O.ID. The energy scale D is 
half the bandwidth. We keep 256 states in the reduced 
density matrix. 



III. COMPARISON OF METHODS 



D. Dynamic density matrix renormalization 



The DMRG introduced by WhiteP^Sl in 1992 is an 
excellent numerical method for one-dim ension al quantum 
systems with open boundary conditionf'^'^ such as the 
SIAM in linear chain representation. The dimension of 
the Hamilton matrix grows exponentially with system 
size. The DMRG provides a well controlled procedure 
to cut off this growth by selecting an optimum basis set 
for the desired states, e.g., the ground state or another 
target state. The optimum basis states are selected from 
the eigenvectors of a reduced density matrix from which 
only the largest eigenvalues are retained. 

We calculate dynamic quantities at zero temperature 
such as the advanced Green functions G{uj—iri) = G^{uj— 
irj) + G^{uj — irj) and G{uj) = lim,,^o+ G{u! — irj) using 
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where AH :— H — Eq. The imaginary part of G< {u — irj) 
provides the spectral densities which we are aiming at. 
Several variants of numerical methods were introduced 
to obtain dynamic quantities, for instance the Lanc- 
zos method^ and the correction vector met hod. 1^211211 
Since the Lan czos method has a limited numerical 
resolution,EiESHlZl for details see the analysis in Ref. [751 
we use the correction vector method which targets the 
ground state |0), the excited state \A) :— A^ \0), and the 
resolvent applied to the excited state. This additional 
targeted state |^±) is called the correction vector 
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Technically, one targets both the real and the imaginary 
part of \£,±). The Green function G'^{uj — irj) is obtained 
from the scalar product 



In this section, we present results for the Hubbard 
model within DMFT using the three impurity solvers in- 
troduced in the previous section. The non-interacting 
density of states (DOS) is given by the semi-ellipse 
PoH = {2/iDn))^l-iLu/Dr. 

The NRG and DMRG calculations were done at zero 
temperature while the ENCA requires a small finite T 
as discussed in Sect. |IIB[ the used values are given for 
each result below. The spectral density p{uj) of the Green 
function is given by p{uj) — ImG(w)/7r and we similarly 
define S(w) := lim,,_j.o+ S(aj — ii]). 



A. Single Particle Dynamics 




G<(a;-iry) = (^|e±) 



(7) 



FIG. 1: (Color online) Spectral function for U — D and /i = 
(left panel) and /i = 0.25D (right panel). While the NRG and 
DMRG results are for T = 0, the ENCA curves are computed 
for T = 0.17D (left panel) and for T = 0.22D (right panel). 
The case ^ = implies half-filling, n = 0.5, while for ^ = 
0.25D the fiUings are riNRG ~ 0.583, tidmrg ~ 0.577, and 
riENCA — 0.573. 

Figure [T] displays a comparison of the local spectral 
densities for a moderate interaction U = D and for 
two different values of the chemical potential. All three 
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FIG. 2: (Color online) Imaginary (top row) and real (bottom 
row) parts of the self-energy for /i = (left column) and 
/i = 0.25D (right column) and U — D. Other parameters are 
as in Fig. [l] 



FIG. 3: (Color online) Imaginary part of the magnetic (top 
row) and charge (bottom row) susceptibility for fi — (left 
column) and /i = 0.25D (right column) and U — D. Other 
parameters are as in Fig. [T] 



methods qualitatively yield the same result with a broad 
central peak and only very small shoulders at energies 
\Ld\ sa D. The latter are precursors of the Hubbard 
bands centered at the energies uj « ±C//2 — /i. The 
NRG and DMRG results agree quantitatively and only 
the shoulders are slightly more pronounced in the DMRG 
curve which is probably due to the lower resolution of 
the NRG at high energies. The central resonance of the 
ENCA curve is narrower and the shoulders of the Hub- 
bard bands are more washed out. The latter feature can 
be attributed to the rather high temperature, T w 0.17D 
for 1.1 — and T « 0.23D for fj, = 0.251?, required to avoid 
the ENCA problems at too low temperatures, see Sect. 
n.B. At large energies, e.g., \uj + /i| > 1.5D, the spectra 
of all three methods agree almost perfectly. 

The corresponding self-energies are shown in Fig. [5] 
For all methods the imaginary part ImS(w) displays a 
quadratic minimum at the Fermi level signalling the va- 
lidity of a low-energy effective Fermi liquid description. 
This implies that the central peak in the spectral function 
of Fig. [T] is essentially due to Fermi liquid quasiparticlc 
excitations. In accord with the Kramers-Kronig (KK) 
relation, Rel](w) behaves linearly at the Fermi level. 
The self-energies of the renormalization group (RG) ap- 
proaches NRG and DMRG agree quantitatively and the 
visible deviations are only due to the different broadening 
procedures used to obtain continuous functions in uj. The 
ENCA approach yields a too steep quadratic and linear 
dependence around the Fermi level in the imaginary and 
real part, respectively. 



B. Collective Modes and Low Energy Scale 

The imaginary part of the local dynamic magnetic and 
charge susceptibilities shown in Fig. |3] shed light on the 



characteristic energies for both types of collective excita- 
tions. The characteristic energy for local charge excita- 
tions is given by the position of the Hubbard bands as 
can be observed in Fig. fl] Consequently Imxchargc('^) 
has a broad peak at roughly uj w 1.1 D. Away from half- 
filling, the asymmetric position of the lower and the up- 
per Hubbard band is reflected by a slightly broadened 
peak. Compared to the charge excitations, the charac- 
teristic energy for local spin excitations is shifted towards 
lower values and the absolute height of Imxmag('^) is 
roughly twice as large as Imxcharge(w). Both these fea- 
tures are signs of the enhanced magnetic Kondo-like cor- 
relations in the system already present for this moderate 
value oiU = D. The overall height of the ENCA suscep- 
tibilities is lower than those of the RG approaches due to 
the nonzero temperature. Additionally, the maxima in 
the ENCA results are shifted to slightly lower energies. 

Increasing the Coulomb repulsion to U = 2D strongly 
enhances the correlations in the system. At half-filling, 
this increase drives the system towards the MIT. In the 
spectral densities depicted in Fig.Hfor half-filling (/i = 0) 
and for finite chemical potential (/z = 0.5D), the Hub- 
bard bands are now well separated from the many-body 
resonance at the Fermi level. The inner edges of the 
Hubbard bands are rather sharp with slight peaks asso- 
ciated presumably with bound trions of a quasiparticle 
and a particle-hole pair.'SSEl] jjj ^j^g ENCA spectra, such 
peaks have also been obtained,'^ but they are not ob- 
served here due to the relatively high temperature. The 
spectra obtained by DMRG and by ENCA are almost in- 
distinguishable at high energies, e.g., for |a;-|-/i| > 1.5D. 
The NRG curve falls off slower for large energies due to 
the limited resolution at large energies mentioned previ- 
ously. 

The deviation between the ENCA and the RG results 
are partly due to the finite temperature to be used in 
the ENCA evaluation. Another part is due to a too low 
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FIG. 4: (Color online) Spectral densities at f/ = 2D and 
chemical potentials /i = (panel (a)) and /i = 0.5D (panel 
(b)). For NRG and DMRG T = holds while for the ENCA 
T = 0.027D ^ O.IST* (/x = 0) and T = O.ID ^ 0.6r* (^ = 
0.5D). For /i — 0.5D, the fillings are nNRG ~ 0.580, nDMRG ~ 
0.592, and wenca = 0.567. 



Fermi liquid scale T* procured by the ENCA. This will 
be shown in the following. 

For further analysis, we extract a common energy scale 
T* from the data of all three methods. In microscopic 
Fermi liquid theory the low energy scale is proportional 
to the inverse mass enhancement!^ determined from the 
quasiparticle weight Z . The energy scale T^ defined in 
this way reads 



T% 



Z D 
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1 - a<^ReI](0) 
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(8) 



and it is shown in Fig. [5] as function of the Coulomb 
repulsion U . For the two values of /i displayed in the 
figure, T^ diminishes with increasing U and vanishes at 
the MIT. For larger values of the chemical potential, e.g., 
1^ ± ^1 > ^, the MIT will not occur anymore due to 
finite doping and T^ will approach a constant value (not 
shown) . 

While the energy scale T^ extracted from the two RG 
approaches agree quantitatively, the ENCA scale follows 
the same trend, but is lower by about a factor of two. 
This in accord with results for the SIAM, i.e., without 
self-consistency, where the ENCA is known to pro duce 
the correct order of magnitude for the Kondo scale.'^SMl 
The ENCA provides the exponential dependence of Tk on 
[/, but the absolute values are slightly too low.'*^ Within 
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FIG. 5: (Color online) Fermi liquid scales r| and T* extracted 
from the quasiparticle weight and from the maximum in the 
spin susceptibility, respectively, as function of U for two values 
of the chemical potential. The magnetic scale T* is rescaled 
by a single factor r| = aT^ = 3.125 • T^ to obtain coinciding 
energies. 



the self-consistency of the DMFT this tendency persists 
and it is slightly amplified. 

Another estimate for the low energy scale can be de- 
fined from the characteristic excitation energies for spin 
fluctuations.^'* We determine a magnetic scale T* from 
the position of the maximum in the local dynamic mag- 



netic susceptibility, T* 



This estimate is equiv- 



alent to the energy scale extracted from the linear slope 
in Imxi„ag(w) for small uj. As can be seen in Fig. p] 
the two scales T| and T* of each method lie on top of 
each other at large U if T* is rescaled by a single factor 
a of order unity. Thus the low energy magnetic exci- 
tations and the single-particle excitations originate from 
the same physical process which is governed by a single 
energy scale. We will call such a behavior 'universal' in 
the context of the present work. For the SIAM such be- 
havior is well known to occur in the Kondo regime. The 
DMFT self-consistency alters only quantitative aspects, 
but no qualitative ones. Hence Fig. [5] indicates universal- 
ity in the metallic phase of the Hubbard model at large U 
where Kondo-correlated quasiparticles dominate the low 
energy physics. 

We observe in Fig. [5|that the magnetic scale T* and the 
single-particle scale T^ differ for small values of U (and 
at large doping, not shown) in analogy to what has been 
found in the SIAM.'^T| and T* differ so that no univer- 
sality can be established. A description of all excitations 
in terms of a single energy scale cannot be maintained. 
The Fermi liquid description is certainly still applicable, 
but all Landau parameters have to be determined inde- 
pendently. 

Henceforth, we write T* — T^ to represent the low 
energy scale and omit subscripts for simplicity. 

The rescaled self-energy T*T,{uj), which determines the 
scattering rate of the Fermi liquid, ^^ is plotted in Fig. [6] 
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FIG. 6: (Color online) Rescaled self-energy T*T,{lo) for U = 
2D and (a) /i = (b) /x = 0.5-D as a function of energy 
measured in units of the low energy scale T* . Within the RG 
approaches we find T* ~ 0.31D and T* « 0.33-D for ^ = 
and /i = 0.5D, respectively. In ENCA we find T* ^ 0.18D 
for and T* ~ 0.2D for jj, — and /i = 0.5D, respectively. The 
temperatures are as in Fig. [4] 



as function of uj/T* for the same parameters as in Fig. HI 
Generally, the agreement between the three methods is 
very good and deviations only occur at large energies, es- 
tablishing the too low energy scale to be the main source 
of discrepancy between the ENCA and RG methods. The 
deviations for oj/T* > 1 observable in panel (b) are on 
the one hand due to th e imp roper description of corre- 
lated valence fluctuationa^EIl and on the other hand due 
to the thermal broadening required in ENCA. 

The rescaled dynamic magnetic susceptibility is de- 
picted in Fig.[7]for two different chemical potentials. The 
peak positions from all three methods coincide but their 
heights differ. For the ENCA this is due to the finite tem- 
perature T oc 0{T*) and the susceptibility is expected to 
increase if T -^ 0. The use of raw NRG data without us- 
ing the equation-of-motion trick^^ might be responsible 
for the discrepancies between the NRG and DMRG sus- 
ceptibilities. 

In the metallic regime, the universality conjectured be- 
fore is supported by the fact that T*ImXrriag(w/T*) ap- 
proaches a universal function for large values of U in the 
half-filled case. This can be observed in Fig.JTJa where the 
susceptibility from NRG is depicted for various values of 
U. For finite /i, however, T*l'niXnia.g{^/T*) continues 
to grow with decreasing T* indicating that valence fluc- 
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FIG. 7: (Color online) (a) Rescaled imaginary part of the 
magnetic susceptibility T*lmx{^/T*) for p = (left panel) 
and jj, = 0.5D (right panel) a,t U — 2D as a function of 
frequency. Parameters are as in Figs. H] (b) T*Imx{'-^/T'') 
vs. Lo/T* calculated with DMFT(NRG) for /i = (left panel) 
and n = 0.25D (right panel) at various values of U. 



tuations modify the low energy physics decisively, thus 
abolishing universality in the lattice model away from 
half-filling. This is in contrast to what is found in the 
impurity model. 

The local charge susceptibility Imxchargc(i^) shown in 
Fig. [8] is strongly suppressed for small energies cj as con- 
sequence of the large Coulomb repulsion. This is partic- 
ularly striking in comparison to the spin susceptibility. 
The characteristic energy scale of the charge susceptibil- 
ity remains set by the interband excitation energy be- 
tween the quasi-particle band and the Hubbard bands 
which is of the order ~ ^ — |/i|. 

In the curves obtained by the RG methods at half- 
filling, peaks emerge at the onset of the interband ex- 
citations as can be seen in Fig. [s] (left panel). They 
originate from the sharp features at the inner Hubbard 
band edges.!^ As in the single particle spectrum, these 
are missing in the ENCA curve due to thermal broad- 
ening. Away from half-filling, the charge susceptibility 
increases due to the increased phase space at low en- 
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FIG. 8: (Color online) Imaginary part of the charge suscep- 
tibility for fi = (left panel) and /i = 0.5D (right panel) at 
U = 2D as function of frequency in units of the half band- 
width D. Parameters are as in Fig. Ill 



ergies. There, the ENCA reveals its limitations in the 
mixed- valence regime since at low temperature its thresh- 
old exponents generate an additional low frequency peak 
for oj — > which is expected to disappear if higher order 
vertex corrections were included. 

For increasing Coulomb repulsion U, the spectral 
weight of ImXchargc is shifted towards larger energies uj 
as illustrated in Fig. 2 in Ref. [T5l In contrast, the po- 
sition of the maximum of the spin susceptibility, which 
sets the scale T* , is shifted towards smaller energies with 
increasing U as shown in Fig. [Tb. As a consequence, the 
energy scales of collective charge and spin excitations are 
clearly separated. This will turn out to be important for 
the observation of kinks. 

Finally, we illustrate the full dependence of the spectral 
densities on the momentum via the bare dispersion e^ 
according to 



P(w,efc) = -Im -—- (9) 

TT cj — 10 — e/j — Ij(oj) 



in the false-color plots in Fig. [9] The separation of single- 
particle energy scales and the influence of particle-hole 
asymmetry can be seen clearly. The almost flat ridge 
around a; = represents the narrow band of heavy quasi- 
particles which is well-separated in energy from the lower 
and upper Hubbard bands below and above uj = 0. The 
coherence of the quasiparticle excitations is lost once ek 
reaches the scale T* where spin fluctuations become im- 
portant. This causes a kink in ReS as will be discussed 
in the following section. At large electron doping shown 
Fig. |9]3, the upper Hubbard band and the quasiparticle 
band merge rendering charge and spin excitations equally 
important for positive energies. However, at negative en- 
ergies, valence fluctuations are suppressed and the sepa- 
ration of energy scales persists. 
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FIG. 9: (Color online) Spectral density p(uj,tk) obtained 
within DMFT(ENCA) as function of frequency uj and bare 
electronic energy e^ for U — 2D at, (a) half-filling {fj, = 0) 
and (b) n = 0.5D. 



IV. KINKS AND COLLECTIVE MODES 

It is well known that the coupling of fermions to ener- 
getically low-lying bosonic modes causes kink-like struc- 
tures in the fermionic dispersion. This picture is based 
on diagrammatic weak fermion-phonon coupling theory, 
see for instance Refs. l86l and [STj and references therein. 
The kink in the fermionic dispersion occurs roughly at 
the bare phonon energy. For stronger coupling, the dia- 
grammatic approach breaks down^^ ^^ and the kink fea- 
ture persists but it does no longer occur at the bare 
phonon energy.M Roughly, the strength of the kink in- 
creases with the coupling between the fermionic and the 
bosonic modes. 

Recently, it was demonstrated that kinks in the elec- 
tronic dispersion are a generic feature of strongly corre- 
lated electron systems without any coupling to external 
bosons.tS Subsequently, it was shown that the kinks oc- 
curring in strongly interacting electron systems can be 
seen as result of the coupling of the fermions to the emer- 
gent collective excitations of magnetic character. Thus 
the system creates its own bosonic modes which in return 
generate the kinks.l^ 

The original argument by Byczuk et al. for the kinks 



was based on the three-peak structure in the spectral den- 
sity p(a;), as shown in Fig. [4] Its essence is as follows: The 
many-body resonance extends around the Fermi level 
a; = from ri_ < to J7+ > 0, where il± are the po- 
sitions of the minima between the many-body resonance 
and the Hubbard bands. 

Numerical results for the positions \fi±\ of these min- 



dressed dispersion 



ima calculated by NRG are shown in Fig. 10 as function 
of the quasiparticle weight Z. The lifting of the degen- 
eracy of |ri±| upon doping is clearly visible. The NRG 
analysis suggests that the power law scaling fl ex Z^/* 
holds for Z ^ (dashed line in Fig. 10 and in its inset). 




FIG. 10: (Color online) Positions of the minima between the 
many-body resonance and the Hubbard bands calculated by 
NRG as function of the quasiparticle weight Z. The inset 
shows the same data in a double-logarithmic plot. The dashed 
line depicts a power law z'^'^^ for comparison. 



Furthermore, the KK relation implies that the real 
part of G{uj) is linear around a; = and that it has 
sharp maxima at positions r± which are of the order of 
the half-width of the many-body resonance. In partic- 
ular \r±\ < \fl±\ holds. For the semi-elliptic DOS the 
self-consistency condition for the effective medium can 
be solved analytically by r(z) = {D/2)'^G{z). For gen- 
eral lattices, using this relation amounts up to approxi- 
mating r[G'(z)] by the linear term of its moment expan- 
sion. Then Eq. (HI) implies kinks in ReS](a;) at w^ with 

|^±| < \n±\. 

An explicit expression for the position of the kinks 
can be obtained by describing the local Green function 
around the Fermi level by 



G(co) 



Z\o 



Lo — loq — i(7 + 7'aj) 



(10) 



The parameters luq, 7, 7', and Zioc are determined 
from the physical quasiparticle weight Z and the non- 
interacting DOS alone. The kink positions wf- are cal- 
culated as the maxima of the second derivative of the 



w^ 



woT 



7 + I'^o 



-7 




(11) 



for details see supplement of Ref. 5, but note the differ- 
ing sign of the inner square root. For the particle-hole 
symmetric case one has 7' = wq = so that ( 11 ) reduces 
to iof = ±(^2- 1)7. For the semi-elliptic DOS j = ZD 
the kinks are thus located at uf = ±(-^2 — 1)ZD. 

Indeed, the kink positions in the numerical data for 
ReS agree nicely with results obtained via Eq. ( 11 ). This 



is demonstrated for /i = 0.25D in Fig. [TT] where the kink 
positions extracted for the DMRG and the NRG impurity 
solver are plotted as function of U/D. Only data for small 
doping and large repulsion is shown because for large 
doping and/or small repulsion the three-peak structure 
of p{uj) is not found so that the above analytical argument 
does not hold and the kinks cannot be resolved. For large 
doping and strong repulsion there still exists one kink in 
ReS for which Eq. ( 11 ) still predicts the correct position 
(not shown). 




-0.2 



NRG I 
NRG ReE 
DMRG ( 
DMRG ReS 



1.8 



2.2 



2.4 

U/D 



2.6 



2.8 



FIG. 11: (Color online) Comparison of the kin k position from 
a numerical analyis of ReE and from Eq. ( 11 1 for /i = 0.25D 



for the two RG impurity solvers as function of U. 

The above analysis due to Byczuk et al. focusses on 
the mathematical structure of the equations defining the 
propagator and the self-energy. In Ref. G\ some of us 
showed for the half-filled case that the characteristic ex- 
citation energy of spin fluctuations agrees with the kink 
energy. Kinks in Rel](w) lead via the KK relation to 
inflexion points in ImS(a;) at the same energies. This 
corresponds to a change in the quasiparticle lifetime 
T ~ l/ImS. Inversely, this implies that humps in ImE 
imply kinks in ReE(cL') via the KK transform. Hence, 
even though no explicit bosonic modes are included in 
the Hubbard model, the emergent collective spin excita- 
tions are responsible for the structures in ImE(a;) and 
thus for the kinks in ReE. 

This argument is also applicable without particle- 
hole symmetry, i.e., for the doped model. The kinks 
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FIG. 12: (Color online) Illustration of the KK transforms of 
two test functions /i and /2, which are depicted in the inset. 
The first, /i displays a kink while the second (/2) does not, 
but has the same average slope for small frequencies. The KK 
transforms are shifted by a^, bi such that KK[fi]{—ai)+bi = 0. 



in Rel](cj) are still associated with additional inflexion 
points in Inil](a;) which are related to changes in the 
relaxation mechanism. To illustrate this view qualita- 
tively, we mimic a kink in Rel](w) by the function /i de- 
picted in Fig. [12] and include for comparison /2 without a 
kink. Then we study the differences in the KK transforms 
which correspond to the imaginary part. While KK[/2] 
is governed by a wide parabola in the range to € (—4, 4), 
KK[/i] displays a noticeable hump starting below the fre- 
quency of the kink. A paraboHc fit would hold only in 
the interval uj € (—1, 1). 

DMFT self-energies computed with the NRG are 
shown in Fig. [13^ . The dashed lines in the upper panel 
indicate the linear fits to ReE(aj) used to determine the 
kink positions. Fig. |13^ displays the same qualitative 
features as Fig. [T2| though they are less pronounced. The 
physical model does not display mathematically sharp 
kinks as the test function /i does. The real part Rel](w) 
in Fig. [13^ displays two kinks. The one at negative fre- 
quencies is fairly clear, the one at positive frequencies is 
fairly weak. Correspondingly, the humps in ImE(a;) are 
clearly visible at negative frequencies, but only weakly 
discernible at positive frequencies. The kinks, which 
mark the beginning of the hurnps, indeed occur at about 
T* = ir*(= iT|), cf. Fig. [sl in agreement with the 
previous finding at half-filling. 

A parabolic description in terms of the Fermi liquid 
scale {ijj/T*Y is possible, but only up to about \ojIT*\ « 
0.3, again in accord with the finding at half-filling.'Sl The 
scattering rate as given by ImE(a;) decreases compared 
to a pure (w/T*)^ behavior with increasing |a;|. 

This picture is consistent with the renormalization 
group flow and the RG fixed points of the effective site for 
a converged metallic DMFT solution. For T* > |cl;| — >• 0, 
the physics is determined by a line of strong coupling 




FIG. 13: (Color online) (a) ReE(cj) (upper panel) and ImE(aj) 
(lower panel) at finite doping as function of loJT* obtained 
from DMFT(NRG). The dashed lines in the upper panel indi- 
cate linear fits used to determine kink positions. In the lower 
panel the parabola (oj/T*)'^ expected from Fermi liquid theory 
is included for comparison. The coherence scale T* ~ 0.16D 
is the low-energy scale for U = 2.6-D, cf. Fig. [5] (b) Same 
as in panel (b), but for [i — 0.5D at f/ = 2D. Note that the 
kink at positive of frequencies is very weak and concomitantly 
ImS(aj) displays parabolic behavior up to oj ~ T* except for 
a very shallow hump. 



(SC) fixed points which des cribes a Fermi-liquid with bro- 
ken particle-hole symmetry.liSESl j^g characteristic energy 
scale is given by T*, and ImS(a;) ex (w/T*)^. With in- 
creasing frequency, however, the system is described by 
the unstable local- moment (LM) fixed point. The dy- 
namic Kondo singlet is broken on a scale Tk oc T* by 
singlet-triplet excitations and the quasiparticles disinte- 
grate at higher excitation energies, leaving a free local 
spin coupled to the conduction band. As a consequence, 
the scattering is reduced and the self-energy is increasing 
much slower than close to the SC fixed point. Spin-flip 
scattering dominates the self-energy in the vicinity of the 
LM fixed point. 

Therefore, the single-particle self-energy retraces the 
crossover from the LM to the SC fixed point. At very 
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high frequencies the magnetic scattering is weak, and the 
physics of the Hubbard model is determined by the local 
charging energies derived from the atomic picture. Hence 
the self-energies depicted in Fig.[6]are low and featureless 
at very high frequencies. At intermediate frequencies, the 
spin-flip scattering provides an additional decay channel 
on top of a weakly correlated, particle-hole asymmetry 
conduction band. This additional decay channel reveals 
itself in l'niE{uj) as the humps at intermediate energies, cf. 
Figs. [12] and [T3| At low energy scales, the spin-flip scat- 
tering is replaced in Iml](w) by the Fermi-liquid parabola 
determined by T* for \oj\ -> 0. 

At larger electron doping, only the lower Hubbard 
band is well separated from the quasiparticle band, see 
Fig. |4j The particle-hole asymmetry and the correlated 
valence fluctuations matter for positive excitation ener- 
gies. The scales for spin and charge excitations are not 
well separated so that both channels contribute to the 
self-energy. The ImS(a;) remains almost quadratic in 
Lu for uj > even on the scale of a; w T* as can be 
seen in Fig. ^ for ^ = 0.51? at C/ = 2D. Only a 
minute hump occurs in the quasiparticle decay rate at 
the spin excitation energy and consequently only a very 
weak kink occurs at positive energies. Even though the 
spin-susceptibility shows a pronounced maximum, the ac- 
cessibility of low-energy charge fluctuations for positive 
energies suppresses the kink, in accord with the two con- 
ditions stated in the Introduction. 

At negative excitation energies, correlated low-energy 
valence fluctuations cannot be excited so that the charge 
energy scale is well separated from the coherence scale 
T*. Thus the kink and the corresponding hump are dis- 
tinct at negative frequencies. 

The above discussion shows that kinks can be directly 
linked to physical processes in the system. As in the sym- 
metric case, the kink positions in the self-energy correlate 
with the positions of the maxima of the spin susceptibil- 
ity, i.e., T* which is shown in Fig. fTllfor two values of 
/x. The values for both quantities from all three meth- 
ods coincide and the small deviations can be understood 
from the strengths and the weaknesses of the methods as 
discussed in Sect. |III A| The kink positions oj^, equal T* 



(a) 0.3 



for r* — >■ 0. This clearly supports the view that the spin 
fluctuations are responsible for the kinks. Deviations oc- 
cur for larger energies corresponding to smaller values of 
U. There, the Fermi liquid description in terms of a sin- 
gle energy scale does not apply anymore since T* and T^ 
stem from combinations of different types of excitations, 
and non-universal valence fluctuations play a role. Addi- 
tionally, kinks are less pronounced for small U and hence 
their positions are harder to determine unambiguously. 
For smaller U the spin fluctuations do not yet behave 
like emergent collective bosonic modes, and the charge 
and spin excitations are not well separated in energy. 

Finally, we study the doping dependence of the kinks 
at fixed interaction. Generic results are depicted in Fig. 
[T5] for hole doping. As observed in Fig. [M] the kink po- 
sitions and energy scale of the spin fluctuations coincide 



3 



NRG 
DMRG 


1 

--X-- 


|m = o 


ENCA 


-■-X--- 




- = t; 












,-*"/ 






JXJ^ 





(b) 



Q 




FIG. 14: (Color online) Kink energies oj^ as function of the 
frequency of the maximum T^/D in the imaginary part of the 
spin susceptibility for (a) /i = and (b) /i = 0.25D. 



uj^ w |w^ I sa T* for small doping 6 £ 0.07. For larger 
doping the particle-hole asymmetry implies that lo^ and 
\ujZ\ differ from each other and hence from T* For hole 



doping, we find ujj > T* > |a;~| but the deviations are 
rather small. Up to an offset, all three energy scales de- 
pend essentially linearly on doping. The energy scales 
rise upon increasing doping. The two panels of Fig. |15| 
compare the doping dependence of the kinks for two dif- 
ferent values of U. Clearly, a larger value of U decreases 
the energy scale of the kinks as one would expect for a 
magnetic energy scale. 

At this point, a comparison to experiment is in order. 
The experimentally best studied strongly correlated sys- 
tems displaying kinks are the superconducting cuprates. 
It is pre sently still debated wheth er th ese kinks are of 
phononicplEll or of magnetic origin.EESl Qur calculation 
based on DMFT and a semi-elliptic DOS is too far away 
from the experimental situation to make a quantitative 
comparison. But it is interesting to note that the kink 
positions observed at low temperatures in underdoped 
high Tf. materials indeed display a linear behavior with 
offset very similar to the one in Fig. 15 Even the num- 



12 



3 



0.25 



0.2 



0.15 



1^ }*'tt ♦■* 




0.05 - 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
5 — I — n 



0.25 



3 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
S=l-n 

FIG. 15: (Color online) Dependence of the kink positions 
cjj^ and of the scale of the magnetic modes on hole doping 
for U — 2D (upper panel) and for U = 2.8D (lower panel) 
obtained from DMFT(NRG). 



bers are in the experimental range^El of about 30 meV 
at zero doping to 120 nieV a,t 6 = 0.15 if we assume 
D = 1.4eV and U = 4eVw 2.8D. The latter number 
is too high by a factor 2 compared to the cupratespbut 
this is not astounding in view of the approximation made 
in the present study. 



We emphasize that we consider only the low-energy 
kinks. These can be expected to be described with an 
effective low-energy single-band Hubbard model which is 
based on the existence of the Zhang- Rice singlets. I^ Any 
statements on the formation of Zhang-Rice singlets and 
on dispersion features at higher energies such as "water- 
falls" which require a multiband description, see for in- 
stance Ref. |95l are beyond the scope of the present work. 
Still, the doping dependence observed here is in accord 
with a magnetic explanation of the kinks in the cuprates. 



V. SUMMARY 

In this paper, we have focussed on the physical origin 
of kinks in the electronic dispersion of strongly correlated 
electronic systems without coupling to external bosons. 
Exemplarily, we have studied the doped Hubbard model 
in dynamic mean-field theory at zero temperature. 

Three different numerical algorithms have been em- 
ployed for solving the self-consistent impurity problem 
for the dynamic mean-field theory and their results have 
been compared. We have found a very good agreement 
between all three impurity solvers. This is quite remark- 
able considering their very different nature and corrobo- 
rates the validity of our evaluations. We distinguish the 
two numerical renornralization group approaches (NRG 
and DMRG) which can be applied directly at zero tem- 
perature and the analytical ENCA which is based on the 
summation of a large subset of diagrams of an expansion 
in the hybridization. 

The NRG is an efficient numerical approach for ar- 
bitrary temperatures and very precise at small frequen- 
cies, but washes out spectral details at higher frequencies 
due to the logarithmic discretization and the concomi- 
tant broadening. The DMRG is the most resource con- 
suming approach because it requires a separate run for 
each spectral frequency. It is used here as a zero tem- 
perature method, though extensions to finite tempera- 
tures are in principle possible. To access exponentially 
small frequencies a logarithmic discretization would be 
necessary.lSfil For equidistant discretization, however, the 
DMRG exhibits a very good resolution at all frequencies. 
The ENCA is fastest in the computation and provides 
data with arbitrary resolution at all frequencies. But it 
suffers from Fermi liquid pathologies below the coherence 
temperature T* which prevent its reliable application at 
very small temperatures. 

With all these three different methods, we have stud- 
ied the kinks in the dispersion of the interacting elec- 
tron system away fronr half-filling. We have established 
that the appearance of the kink is linked to dominant 
spin fluctuations at low energies. In particular, the po- 
sition of the kinks in energy is intimately linked to the 
spin-fluctuation scale T* where the magnetic suscepti- 
bility exhibits its maximum. Additionally, the occur- 
rence of kinks requires a substantial energy separation 
between the charge fluctuation scale U/2 — |/i| and T* 
i.e. f7/2— |/x| 3> T* . This is clearly the case in the strongly 
correlated regime at large values of [/ ^ 2D = W where 
W is the bandwidth. Therefore, the low-lying bosonic 
modes of the electronic system modify the electronic dis- 
persion on the scale T* at sufficiently strong electron- 
electron interaction U. 

Universality, i.e., the possibility to describe the low- 
energy dynamics in units of a single energy scale, is only 
observed close to half-filling where the scale is set by 
the spin-fiuctuation scale T^. For large doping and/or 
weaker interaction, there is no clear separation between 
the spin- and charge-fluctuation scale and universality is 
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lost. 

Finally, we have compared the doping dependence of 
the kinks in the single-band Hubbard model with kinks 
measured with angle-resolved photoemission in planar 
cuprates. In spite of the much higher complexity of the 
cuprates compared to our model study, the energy scales 
are qualitatively reproduced: Magnetically induced kinks 
evolve from 30 nieV to about 120 meV with increasing 
hole doping for U = lAW. 

In conclusion, we extended the view that emergent 
modes of the electronic system can generate kinks in 
the electronic dispersion in analogy to materials with 
strong electron-phonon coupling even away from half- 
filling, provided (i) that there is a significant energy sepa- 
ration between the high and low energy scale, T*/D <C 1 
andwhich implies that we are in the strong coupling 
regime U/W ^ 1, and (ii) that the charge energy scale is 



much larger than T* Note that in the studied strongly 
correlated system the ratio T*/D takes the role of the 
ratio of the Debye frequency over the Fermi energy in 
coupled electron-phonon systems. We expect that this 
a generic feature and very similar results will apply in 
many other related systems with clearly separated en- 
ergy scales and strong emergent modes. 
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